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ABSTRACT 



A review of the current theory and methods for the computa- 
tion of maximum flow in networks is presented along with a simplified 
method for determination of a feasible flow in networks with upper 
and lower bounded arcs. A computational procedure is presented which 
is used to calculate the maximum flow for a general network. The 
network is reduced to an equivalent basic network. An associated 
network is used to compute a feasible, then the maximum flow for 
the basic network. A computer program is included for use in 
computation of maximal flows in large networks. 
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CHAPTER I 



INTRODUCTION 

The growing uses for mathematical programming models In trans- 
portation and communication networks have prompted study in the devel- 
opment of computational methods to determine maximal network flows. 

The subject of this thesis is the mathematical problem of determining 
maximal steady state flows in networks which are subject to capacity 
limitations on the arcs and nodes of the network. 

An efficient and widely used method to compute maximal flow 
was developed by Ford and Fulkerson [7]. This method along with a 
review of the current theory will be presented in Chapter II in 
order to give the reader the theoretical foundation upon which the 
remainder of the thesis is based. Chapter III introduces a new 
method which utilizes the Ford and Fulkerson algorithm to compute the 
maximal flow in networks with non-zero lower bound limitations on the 
flow in the arcs. Since many networks of interest are more general 
In nature, the generalized network is introduced in Chapter IV along 
with a computational procedure for the determination of maximal flow. 
To illustrate the techniques proposed, an example is included. The 
procedure can be used to solve small problems by hand, but the amount 
of work increases rapidly with the size of the network. Therefore, a 
computer program which may be used to solve larger problems is 
included in the appendices. A brief description of the program is 
found in Chapter IV. 
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Definitions and Symbols 

In order to establish a common understanding as to the exact 
meaning of various terms and symbols used to present the material of 
this thesis, an initial set of definitions is presented. Other terms 
will be defined throughout the material as needed to facilitate the 
presentation. 

Define a network , G, as a graph which can be represented on 
a plane in such a way that the set, N, of nodes are distinct points 
and the set, A, of arcs are simple curves which connect two distinct 
nodes. Furthermore, no two arcs can meet except at the nodes which 
are their extremities. The graph contains no loops. 

Define io"(N) to be the set of arcs or flows incident to a set 
of nodes, N, and oo + (N) to be the set of arcs or flows incident from N. 

A directed arc, a 4 _., Is defined as an arc incident from a node 
i and incident to a node j. 

An undirected arc connects node 1 and node j without having an 
orientation. 

Flow can be thought of as a value of the steady state rate of 
movement of a homogeneous commodity along a path or channel. Positive 
flow in a directed arc, a. . , will move with the orientation of the 

* »J 

arc and will be denoted by the symbol, y^ y Positive flow In an 
undirected arc, denoted y^ y may move in either direction but not 
in opposite directions at the same time. 

To each arc, a 4 . £ A, there will be associated two real 

» *J 

numbers, b^^. , and c^ j, which represent, respectfully, the minimum 
and maximum allowable flow In that arc. The bounds associated with 
undirected arcs, a^ y are denoted b^ j, and c^ j. 
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An arc Is saturated If j * j. 

Define as the set of flows in a network G. 

Initially, the capacities of the nodes will be assumed to be 
Infinite. The nodes will be classified according to the value of the 
difference of the sum of the flows Into the node minus the value of 
the sum of the flows out of the node. A node Is designated as a 
source node If this difference Is negative, and as a sink node If this 
difference Is positive. For all Intermediary nodes, the sum of the 
flows into a node must equal the sum of the flows out. The set of 
source. Intermediary and sink nodes will be denoted as S, I, and T 
respectfully. 

Define F as the value of the flow In a network. 



F = Z y, t 

< T > 

A flow Is feasible If and only If: 

b 1,j - y 1J - c 1,j V a 1,j - A 

(2) F - 2 + y , 

^ < s > ,;l 



1 - y 1 t " F 

y 1 , t - w < T > ’ 



= 0 



(3) 



2 y " 
y^jCufO) y 1 f jCuT(j) 






Ucl 



Define a chain as a sequence of arcs such that each arc, a^ j 
1,j c I, Is connected to an arc a^ ^ at node 1 and to arc aj -j at node j. 
A cycle Is a sequence of arcs where: 

(1) The sequence Is a chain. 

(2) The sequence does not use the same arc twice. 

(3) The Initial and terminal nodes of the chain coincide. 
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A cut Is defined as any set of directed arcs containing at 
least one arc from every chain of positive capacity joining the 
source to the sink. The value of the cut Is equal to the sum of the 
capacities of the arcs of the cuts. 

Basic Network 

Much of the Initial development of the material In this thesis 
deals with computation of flows In simple or basic networks. There- 
fore, for the purpose of this paper, a basic network, G, Is defined 
as a network which has the following properties. 

1. The network has a return arc, a n ^ , which Is added to the 
network for computational purposes only. 

2. Node 1 has only one arc, a* Incident to It. Node 1 Is 
designated as the source. In the basic network there Is 
only one source node. 

3. Node n has only one arc, a n ^, Incident from It. Node n 
Is designated as the sink. In the basic network there Is 
only one sink node. 

4. For each arc, a^ j, the associated c 1 j, Is positive and 

b i.r°- 

5. Two distinct nodes (1,j) may be connected by only one 
directed arc. The network contains no undirected arcs. 
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CHAPTER II 



REVIEW OF THE LITERATURE 

In order to acquaint the reader with the subject of maximal 
flow in networks, a review of the current theory and primary combina- 
tional techniques is presented in this chapter. 



Basic Network Theory 

A great deal of theoretical work in networks has already been 
done. A few of the relevant theorems are presented here in order to 
provide justification and a rationale for the combinatorial techniques 
and algorithms that are presented in later sections. 

Theorem 1. Opposite directed flows on the same arc cancel [4]. 
This theorem means that given an undirected arc, a. ., with flows 

■ >J 

y. . 0 and y. . >.0, these two flows can be replaced by the flows 

y 'i,j andy 'j,i 

y 'i.o w h,i - min( *i,j * yj.i> 

Theorem 2. A set Y with F>0, satisfying the capacity con- 
straints and node conservation equations can be decomposed into a set 
of positive chain flows from the source to the sink and a set of 
circular flows such that the direction of positive flows in any 
common arc is the same for all chains [4], 

Theorem 3. There exists a positive flow from the source to 
the sink if there exists at least one chain of arcs with positive 
capacity from the source to the sink [4], 

Theorem 4. A flow F° is maximal if and only if the maximal 
flow is zero in a second network formed by replacing c. . by c. . - 

*1 1 jJ 

y i >j v a i .£ A - 



n 



Theorem 5. Given any partition of nodes into two sets, where 
the first set includes the source node and the second set the sink, 
then a feasible solution F is maximal, if every arc, a. ., that joins 

■ sJ 

a node in the first set to a node of the second is saturated [4]. 

Theorem 6. The maximum flow in a network is equal to the 
minimum cut. 



Tree Method 

Dantzig [4] has developed a systematic procedure for com- 
puting the maximal flow in a network by using chains connecting the 
source and the sink of the network. A knowledge of the tree method 
will give a better understanding of the Ford and Fulkerson algorithm 
which is presented in the next section because the two methods are 
identical in principle. 

A tree [1] is defined as a connected graph with at least two 
nodes and no cycles. 

Consider a basic network such that: 

b i,j = 0 • c i.j ^ 0 Va i,j^ A 



The arcs may be directed or undirected. Initially all arcs 
are unmarked and the feasible flow is zero. In order to determine 
the maximum flow: 

1. Create a tree such that the arcs are unsaturated. 

2. Select two sub trees - one , branching out from the source, 
node 1, and the other, T n , branching out from the sink, n, such that 
each node is reached by just one arc of the tree. 

3. Choose one unsaturated arc which connects the two trees. Thus, 
there will be just one chain from node 1 to node n. If none can be 
found, go to step 7. 



12 



4. Determine the amount that the flow can be increased along the 
chain from node 1 to node n, by the following procedure. Determine 
the amount that the flow can be increased, AX, where AX = min (c. . - 
y. Yr, i)* Flow along the chain directed toward the sink is 

1 9 J K J I 

designated y. . and flow directed toward the source along the chain 
■ >J 

is designated y^ -j. There are two cases: 

(a) AX > 0, increase the flow by AX along the chain, go to 
next step. 

(b) AX = 0. This means one arc in the chain is saturated 
by a flow directed toward the sink, go to next step. 

5. Mark one of the saturated arcs. 

6. Eliminate the saturated arc from the chain. Thus, there are two 
sub-trees. Return to step 3. 

7. Check all marked arcs joining T-j to T . If all these arcs have 
flow from T-j to T n then the optimum has been reached. If any arc 
has flow from T n to T^ , use this arc to connect T-j to T n and proceed 
with step 4. 

Ford and Fulkerson Algorithm 

Given the basic network, the most direct method developed to 
date of computing the maximum flow is the algorithm developed by Ford 
and Fulkerson. This method is the one most cited in the literature 
concerning maximal network flow. The procedure is outlined below. 

1. Begin with a feasible flow in 0 in all arcs. 

2. Label each node as follows: 

a. Label node 1, with the number 1. 
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b. If node i is labeled and node j is not labeled, label 

node j with the number (i) if: 

(1) a. . exists and y. . < c. , 

1 5 J I 1 sJ 

(2) a. . exists and y. . > 0 

J 9 ' J 9 ■ 

c. If the sink, node n, becomes labeled by this procedure, 
then the flow from node 1 to node n can be increased. 

3. If node n is labeled, construct a simple chain from node 1 
to node n by back tracking from n to 1 using the labels 

on the nodes. 

a. The chain created will have the property that all arcs, 
a. .on the chain directed toward the sink have a flow 
0 < y. . < c. . and all arcs a, . directed toward the 
source will have a flow y^ ^ > 0. 

b. The flow can be increased by an amount AX where: 

AX = rnintc, j - y,,j. y,, k > * ,k in cha1n - 

4. Increase the flow along the chain by an amount AX and 

return to step 2. 

5. This procedure is repeated until the sink cannot be 
labeled. Flow is maximized when this is the case. 

Now consider the more general problem in which the lower bounds 

on the directed arcs are no longer assumed to be zero. Thus the flow 

in each arc is bounded both above and below. 

Such that: - “ b. . <_ y. . £ c. . £ + °° a. . c_ A 

Theorem 7 r In the case of networks with lower bounds on the 
arcs, C is defined as the subset of N such that if node 1 ^ C the node 
n ^ C or if node 1 £ C then node n £ C. 
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A 

1 . 

2 . 



flow y. . exists and is feasible if and only if: 

1 >J 

b U- y U- c i,o Va i,j^ A 

E c. . > E , b . . 

a, , c of (C) 1,J a, , c w (C) 1 ,J 



i >J - 



i ,J - 



This more general type of problem can be satisfactorily worked 
by using the Ford and Fulkerson algorithm if a feasible flow can be 
found. 



The Ford and Fulkerson algorithm can be used to solve for a 
maximum flow in an associated network G' which is derived from the 
original network, G. This method is outlined below. 

1. Create an associated network G' such that G 1 contains 
the nodes and arcs of the original network G with arc 
capacities c'. . = c. . - b. . . The lower bound, 

• >J • >J I >J 

b'. . equals zero. 

■ >J 

2. Add a new source, labeled 0, and a new sink, labeled m, 

(m = n + 1 ) , to G' . 

a. If the arc, a. ., exists, then construct two new 

* >J 

arcs according to the value of the lower bound, b. .. 

* >J 

If b. . > 0, then construct the arcs, a'. m of capacity 
I »j i ,m 

c'. = b. . ana a* . with capacity c' . = b. . : 

i .m i,j o,j K J o,j i , j 

If b. . < 0, then construct the arcs, a' .'’V/ith 
i »J 0,1 

capacity c' . - - b. . and a ' . with capacity 

o > i i > j J 

c' . = - b. . . 

J >m i ,j 

b. The Ford and Fulkerson algorithm can then be used to 
determine the maximum flow in G' from node 0 to node m. 

3. Two cases can result from step 2b. 

a. Case one: If the flow in each arc, a' . and a 1 . 

o,0 J»m 
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j c N 1 , equals the capacity of each arc, then the flow 
is feasible for the network G and the problem has a 
solution. 

b. Case two: If for any arc, a' . or a' . j £ A' , 

o j j J 

the flow in that arc does not equal the capacity of 
the arc, the problem has no solution. 

4. If a feasible flow for G can be found by the above pro- 
cedure, transfer the flow from G' to G by using the 
following transformation: 



y. . = y' . . + b. 
^i.J y i.J i,J 



V a. . c A 
1 »J 



5. Ford and Fulkerson's algorithm is then used to determine 
the maximum flow using the feasible flow as a starting 
basis. 
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CHAPTER III 



A SIMPLIFIED METHOD FOR COMPUTATION OF MAXIMAL FLOW 
To compute a maximum flow in networks with non-zero lower 
bounded arcs, an initial feasible flow must be determined. The method 
presented in the last chapter produces a complex associated network, 
G 1 , even with a relatively simple network, G. The principle advantage 
of the method presented here [10] is that an initial feasible flow 
can be determined by use of a more simple network, G 1 , than in the 
previous method. 



Proposed Algorithm 

The development of the algorithm is based upon a consequence 
of the flow balance equations. In order to provide a rationale for 
the procedures presented, the development of the associated network 
will be presented now. 

Given a basic network G where the flow in each arc is bounded 
as follows: 



- * < b . . < y. .<c. .+oo 
- i ,J - J i >J - i,J 



a i ,j £ A 



The problem of determining the maximum flow in G can be stated as 
follows: 



Maximize: 
Subject to: 



Vi 

s y 



. , - X , y, , = 0 V i c N 
J ’ y. . cuj (i) 



y. . £ of (i) J> y, ,cu 

J > 1 1 > J 



Va i,j £A 

Now transform the entire network G into an associated network 

G' using the change of variables: 

y. . = y 1 . .+b. . Va. . c A 
*1 ,J y i »J i >J i ,J - 
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